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ABSTRACT 

There are several different methods to determine the individual supernovae (SNe) initial 
explosion energy, here we derive the average or typical explosion energy of shell-type supernova 
remnants (SNRs) in a particular way. By solving a group of equations pertaining to shell-type 
SNRs at the same stage we obtained some physical parameters, e.g. the distance (d), evolved 
age (t), etc.. Assuming series of different SN initial explosion energies ranging from 10 48 ergs 
to 10 53 ergs, we derived series of distance and age parameters with which compared already 
known ones. Thus the most likely value of the SNe initial explosion energy is obtained when the 
deviation is least, which equals to about 10 51 ergs, in good agreement with the undertook value. 



Subject headings: supernovae: general — energy — value 



1. Introduction 

We have already known the initial kinetic en- 
ergy (Eq) of lots of Galactic SNRs through vari- 
ous estimate methods (tabled]). These SNRs em- 
body shell-type and composite-type ones. Fig. Q] 
denotes the number distribution of 44 such rem- 
nants. It seems that they concentrate to about 
10 50 ~ 10 51 ergs. Therefore we can significantly 
take the average initial energy (Eq) as a typical 
SNR physical parameter which is also able to be 
determined by our unique means. 

Our mathematical method to derive SNe initial 
energy is rather different from others. One first 
method of others, for example, by E = \M su v 2 , 
M su is the swept- up mass of remnant shell expand- 
ing into interstellar media (ISM) , v is the velocity 
of shock wave of remnant, the initial explosion en- 
ergy of SNR G180.0-1.7 is thus obtained (Braun 
et al. 1989). Sun et al. (1999) calculating their 
detected ASCA data plus ROSAT data has de- 
rived the initial energy of G327. 1 — 1.1 as the fit- 
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ting result. After knowing the SNR G299. 2-2.9 
radius (R) value, the particle density (no) and 
the age (t), Slane et al. (1996) get Eq value by 
Eq w 340i? 5 n ^ 2 xl0 51 ergs. Bambaet al. (2001) 
have got the Eq value by assuming a thin thermal 
NEI plasma model plus standard Sedov model. 
And so on. All of them obtain the initial energy 
of a individual SNR through physical means. Here 
in the paper we get SNe average initial explosion 
energy by statistical method. 

Galactic supernova remnants are classified 
into three types: shell-type, Plerion-type and 
composite-type. Our work here merely include the 
shell-type remnants. Moreover, shell type SNRs 
usually have four evolution stages: the free ex- 
pansion phase, the Sedov or adiabatic phase, the 
radiative or snowplough phase and the dissipation 
phase. Nearly all the observed SNRs are in the 
adiabatic phase, or in the 3rd. And almost none 
is detected in the 1st and 4th phases. 

In the paper, numerical analysis for the most 
likely value of the SN initial energy is described 
in section [21 of which they are made separately at 
adiabatic-phase and radiative-phase for shell-type 
SNRs, and both made comparison by distances 
and ages. In the last section we discuss and 
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Table 1 

List of the initial explosion energy E of 44 Galactic shell-type or composite- type SNRs 

OF WHICH THEIR VALUE WAS SOMEWHAT WELL DETERMINED BY VARIOUS METHODS. 



Source 


E x 10 51 


Ref. 


Source 




E x 10 51 


Ref. 


— 


ergs 


— 


— 






ergs 


— 


GO.0+0.0 


>40 a 


KF96 


G180 


.0- 


■1.7 


0.24 


BGL89 


G4.5+6.8 


0.4/0.5 b 


BKV05 


G184 


.6- 


■5.8 


0.015 


GAT04a 


G18.8+0.3 


0.01 


D99 


G261 


.94 


-5.5 


0.29 


CD80 


G18.9-1.1 


0.08-0.18 


H04 


G263 


.9- 


3.3 


1-2 


GAT04a 


G24.7+0.6 


0.066 


L89 


G290 


.1- 


■0.8 


0.8 


S02 


G28.6-0.1 


0.9 


BUK01 


G291 


.0- 


■0.1 


0.25? 


HHS98 


G29. 7-0.3 


2 


HCG03 


G292 


.0+1.8 


0.18 


GAT04b 


G31.9+0.0 


0.3-1.4 c 


C05 


G292 


.2- 


0.5 


1.0 


GAT04b 


G39. 7-2.0 


1 


MS96 


G296 


.1- 


0.5 


0.23 


HM94 


G41. 1-0.3 


1.2 


SDP05 


G296 


.5+10.0 


0.2-0.6 


GAT04b 


G74.0-8.5 


0.24 


BGL89 


G299 


.2- 


2.9 


0.12 


SVH96 


G82.2+5.3 


0.17 


M04 


G312 


.4- 


0.4 


0.6 


CB99 


G93.3+6.9 


0.39 


LRR99 


G315 


.4- 


■2.3 


0.66 


GAT04b 


G94.0+1.0 


>0.27 


F05 


G320 


.4- 


■1.2 


1-2 


GAT04b 


G106.3+2.7 


>0.07 


KUP01 


G326 


.3- 


■1.8 


1.0 


GAT04b 


G109. 1-1.0 


1-10 


GAT04a 


G327 


.1- 


■1.1 


0.23 


SWC99 


Gill. 7-2.1 


2-3 


V06 


G327 


.6+14.6 


1.0 


GAT04b 


G116.9+0.2 


0.1 


CHP97 


G347 


.3- 


0.5 


1.0 


MTT05 


G119.5+10.2 


0.03 


GAT04a 


G349 


.74 


0.2 


0.5 


L05 


G120. 1+1.4 


1.16 


WHB05 


G352 


.7- 


■0.1 


0.2 


K98a 


G126. 2+1.6 


>7 


B05 


G357 


.74 


■0.3 


0.11 


L89 


G132. 7+1.3 


0.31 


GAT04a 


G359 


.0- 


0.9 


0.21 


L89 



a Notes: In this case we just make use of the minimum value for our purpose. 
b Notes: In the case we take the average value. 
c Notes: The same as above. 
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summarize our results. 

2. Numerical Analysis 
2.1. At adiabatic phase 

Let us list the group of equations for shell-type 
supernova remnants at the second stage as fol- 
low (Wang & Seward 1984, Koyama & Meguro-Ku 
1987, Bignami & Caraveo 1988, Xu et al. 2005), 

D pc = 4.3 x 10- 11 ( ^ ) 1/5 ^/5 (1) 



E(D) - 1.505 x l - 19 -p^ 

arcmin 

= 2.88 x 10- 14 £r 3 - 8 nl- 3 ( 2 ) 

( lO^ergs ) = 5.3 X 10 'll^-s^l s- 1 

/ n \ 312 

x ( % ) (3) 

Here, D pc is the SNR diameter in unit of pc, t yr 
is the remnant age in year, n is the ISM electron 
density in cto~ 3 , S\ghz is the detected fluxes of an 
SNR in Jy at 1 GHz, 6 arcmin is the observational 
angle in arcmin, v — %P is the velocity of shock 
waves in km s^ 1 . And we know tan ( 6arc ™ in ) = 

2^") where, d pc is the distance to a remnant in 
pc! 

From these parameters above, the fluxes Sighz 
and observational angle 9 arcm i n is the detected 
value for each SNR (table [2]) of which we regard 
these remnants evolving at the Sedov-phase since 
their diameter less than 36 pc. And the SNRs ex- 
plosion energy E — £ x E is (£ = 1, 2, 3, 10 5 ) 
is an assumed value. But the diameter D pc (and 
distance d pc ), age t yr , velocity Wfe ms -i and elec- 
tron density n cm -3 is unknown and to be derived. 
Parameters of the distance d pc = d our (^i), and 
the age t yr = t our (£,i) (£ = 1, 2, 3, 10 5 , i = 1, 
2, 3, 37) will be adopted in our paper. But we 
do nothing with v and n. 

For a certain supernova remnant (i) (i = 
1, 2, 3, ...,37) (table [2]) and assumed series of SNe 
initial explosion energies Eq = £ x i?4s(^ = 



1, 2, 3, 10 5 ) in the unit of 10 48 ergs, we can ob- 
tain the remnant distance d our (^,i)(= d pc ) and 
the age t our (£,,i)(= t yr ) values by solving the 
group of equations above. Then we compare them 
with the already known parameters d trU e(i) and 
ttrue{i) (i — 1, 2, 3, 37) listed in table [2] in or- 
der to derive the most likely original energy of 
remnants. Thus the explosion energy was derived 
when the deviation in comparison is least. 

The group of equations are not strictly correct 
as not to be figured out mathematically, but they 
are correct enough for us to confirm the SNe initial 
energy (E ). 

Figure [2] shows that when Sigh z,il 'Sigh z,i = 
9li '61 = 1 for both SNR1 and SNR2, here 0?(i = 
1,2) is the visual area of the remnant, their ra- 
dio surface brightness (S) can be the same to 
each other. Thus the remnant diameter (D) and 
distance (d) value will be uncertain according to 
equation (2). But fortunately the true reality 
will never take on this case because one can see 
E(D) ~ £>~ c 3 8 from formula (2) and not ~ D p ?. 
Therefore we are able to uniquely determine the 
SNe initial kinetic energy (Eq). 

Many of the radio SNRs have more than one 
published value for distance and age in table [2] 
and [3] For these, we either chose the most re- 
cent estimates or used an average of the available 
estimates, or the most commonly adopted value. 

2.1.1. Comparison by distances 

We can compare these resolved distances 
{d our (^,i)) above with the already known ones 
(dtrue(i)) listed in table [2] by 

37 
i=l 

n 
i=l 

(£=1,2,3,...,10 5 ) (4) 

For example, one gets $(10, d) — 0.126, when 
£ = 10, and $(100, d) = 0.184, when £ = 100. 

Figure [3] shows that the most likely value of 
supernova initial explosion energy (Eq) derived by 
this method for the shell-type remnants at Sedov- 
phase equals nearly to 0.23 x 10 50 ergs. 

There are some different methods to derive 
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Table 2 

List of the distance (d), Age (t) and some other physical parameters of 37 shell- type 
Galactic SNRs of which their diameter is less than 36 pc. 



Source 


l-true[l) 


Utrue\ 1 ') 
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pc 
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99 
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Q 
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4: 


f; 
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PS9 PnzLn 


p^i Q-uo n 




79nn 




7-v- c i 

( AO 


94 




uOZ.o — U. 1 




71 nn 

1 1UU 
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1 7 


1 1 
1 1 


ivyo u 


poo fi_un 1 

vjOO.Otu. 1 


yuuu 


7snn 

/ OUU 


9"? 


i n 


99 
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1UUU 
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1 1UUU 


99 
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Q 
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QQK PnzLa 
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1UUUU 


1 ^nn 

lOUU 


Q 
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99? 
zz . 
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y 




P7zl & ^ 




zinn 


93 
zo 


93nvi fin 


91 n 

Z 1 u 




P7S 9_l_9 1 


^nnnn 


1 ^nn 

lOUU 


9fi 
zu 


fin 


^zin 

o^u 




P&zl 9 OR 


1 1 nnn 

1 1UUU 


zi^nn 


93 
zo 


9n-iri fi 


1 1 


M^JSn M77 PnzLci 


PSQ 0-L/L 7 


i onnn 
ryuuu 


snn 

OUU 


94 
Z^ 


i 9n^on 
izuxyu 


99n 
zzu 


T A Qfi 


vjyo.o-ro.y 


^nnn 

OUUU 


99nn 
zzuu 


10 


97v9n 
Z ( xzu 


Q 

y 


T QQ r*n4a 


PQQ 7 9 
uyo. / — u.z 




1 ^nn 

lOUU 


oo 


sn 

oU 


fi^ 
oo 


U 1VOUZ 


p i OQ 1 1 (1 


1 7nnn 
1 / uuu 


^tnnn 

OUUU 


94 
Z^l 


98 
Zo 


9n 

ZU 


VTTQ^ TTTTirSI PnzLa 
rxiyo, xixivoi, vju^d 


P 1 1 1 7 9 1 




3/Lnn 

OIUU 


o 





979n 
z I zu 


ir vui 


P 1 1 zl ^J_0 ^ 


zii nnn 

^rlUUU 


7nn 

f uu 


IO 




R 
U 


MRPn9 PnzLa 


p 1 1 f\ i 


oennnn 

ZoUUUU 


i finn 

lOUU 


39 

oz 
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OUXOU 


1 1 
1 1 


Xl-Dol, vjrU^td 


pi i a q j_ n 9 


zizinnn 


1 finn 

1UUU 


1 fi 


o^t 


Q 

y 




G120.1+1.4 


410 


2300 


5 


8 


56 


H90, G04d 


G260.4-3.4 


3400 


2200 


35 


60x50 


130 


B94, RG81 


G272.2-3.2 


6000 


1800 


8 


15? 


0.4 


D97 


G284.3-1.8 


10000 


2900 


20 


24? 


11 


RM86 


G299.2-2.9 


5000 


500 


2 


18x11 


0.5 


SVH96 


G309.2-0.6 


2500 


4000 


16 


15x12 


7 


RHS01 


G315.4-2.3 


2000 


2300 


28 


42 


49 


DSM01, G04a 


G327.4+0.4 




4800 


29 


21 


30 


SKR96, WS88, G04a 


G327.6+14.6 


980 


2200 


19 


30 


19 


SBD84, G04a 


G332.4-0.4 


2000 


3100 


9 


10 


28 


CDB97, MA86, G04a 


G337.2-0.7 


3250 


15000 


26 


6 


2 


RHS01 


G337.8-0.1 




12300 


27 


9x6 


18 


K98b 


G346.6-0.2 




8200 


19 


8 


8 


K98b, D93 


G349. 7+0.2 


14000 


14800 


9 


2.5x2 


20 


RM01, G04a 


G352. 7-0.1 


2200 


8500 


17 


8x6 


4 


K98a 



a Notes: Diameters were calculated by using the distances together with the angular 
sizes in Green (2006) catalogue. 
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Fig. 1. — the number distribution of Galactic su- 
pernova remnants corresponding to their progeni- 
tor initial kinetic energy (Eq). 




observer 



Fig. 2. — The plot shows for both SNRs their ra- 
dio surface brightness (E) can be equal to each 
other when S 1GHz ,2/ 'S\ghz,x = = 1- Here 

9f(i = 1,2) is the visual area of the remnant. Thus 
we are not able to uniquely derive the SNR diam- 
eter (D) and distance (d) value. But fortunately 
this case will never occur since E(D) ~ D~^ 8 and 
not ~ -Dj7 c 2 - Therefore we can obtain the just rem- 
nant distance (d) by solving the group of equations 
offerred in text. 



the distance value (d) of supernova remnants (Xu 
et al. 2005) in Galaxy, and majority of these 
distances are not taken advantage of the (Eq) 
value obtained by other authors before. These al- 
ready known distances were derived by the opti- 
cal proper motion/ velocity, HI absorption, associ- 
ation with CO/HI/HII region, HI column density, 
pulsar parallax and optical velocity, or by various 
methods (e.g., Green 2004; Xu et al. 2005), but 
not by T,-D relation or Sedov solutions. Therefore 
through comparison by distances to confirm SNe 
explosion energy (Eq) is to some extent reasonable 
and causes almost no self-contradiction. 

2.1.2. Comparison by ages 

Similarly we can compare these resolved ages 
(iow(£j *)) above with the already known ones 
{ttrueij)) in table [2] by 

37 

= Z)(W(f,0-ttr«e(*)) 2 
i=l 

n 
»=1 

(£ = 1,2,3,...,10 5 ) (5) 

For example, one gets $(10, t) — 0.532, when 
£ = 10, and $(100, t) = 0.507, when £ = 100. 

Figure [4] shows that the most likely value of 
supernova initial explosion energy (Eq) derived by 
this method for the shell-type remnants at Sedov- 
phase equals nearly to 7.0 x 10 50 ergs. 

There are many measures to obtain the age 
values (t) of remnants (Xu et al. 2005), and a 
majority of these ages are not made use of the 
Eq value derived before. For example, if a rem- 
nant is associated with a pulsar, we can estimate 
its age by using the neutron star characteristic 
age obtained from the rotation period of the pul- 
sar (P) and the rate of change of period (P) by 
t = P/2P (Gotthelf et al. 2000). For SNRs 
with a known radius (R) and thermal tempera- 
ture (T) taken from X-ray data, one can obtain 
the age by t = 3.8 x 10 2 P pc (fcT),; c \/ 2 yr (Seward 
et al. 1995). We can also calculate the SNR age 
byf w 40000P~ 1 ' 5 ^ b ~ ' 5 yr, when a remnant has 
its spectrum showing the usual break at frequency 
Ub due to synchrotron losses in a magnetic field B 
(Bock et al. 2001). Therefore through compar- 
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Table 3 

List of the distance (d), Age (t) and some other physical parameters of 20 shell- type 
Galactic SNRs of which their diameter is larger than 36 pc. 



Source 




dtrue (0 


Dia. a 


size(#) 


SlGHz 


Ref. 




yr 


pc 


pc 


arcmin 


Jy 




G8. 7-0.1 


15800 


3900 


51 


45 


80 


G96 


G18. 8+0.3 


16000 


14000 


57 


17x11 


33 


D99,G04a 


G49. 2-0.7 


30000 


6000 


52 


30 


160 


KKS95, G04a 


G55. 0+0.3 


1100000 


14000 


71 


20x15? 


0.5 


MWT98 


G65.3+5.7 


14000 


1000 


78 


310x240 


52 


LRH80, R81 


G119. 5+10.2 


24500 


1400 


37 


90? 


36 


MOO 


G127.1+0.5 


85000 


5250 


69 


45 


13 


FRS84 


G132.7+1.3 


21000 


2200 


51 


80 


45 


GTG80, G04a 


G156.2+5.7 


26000 


2000 


64 


110 


5 


RFA92 


G160.9+2.6 


7700 


1000 


38 


140x120 


110 


LA95 


G166. 0+4.3 


81000 


4500 


57 


55x35 


7 


L89, KH91, G04a 


G166.2+2.5 


150000 


8000 


186 


90x70 


11 


RLV86 


G182.4+4.3 


3800 


3000 


44 


50 


1.2 


KFR98 


G205. 5+0.5 


50000 


1600 


102 


220 


160 


CB99 


G206.9+2.3 


60000 


7000 


102 


60x40 


6 


L86 


G266.2-1.2 


680 


1500 


52 


120 


50 


K02, AIS99 


G296. 5+10.0 


20000 


2000 


44 


90x65 


48 


MLT88 


G296.8-0.3 


1600000 


9600 


47 


20x14 


9 


GJ95, G04a 


G321.9-0.3 


200000 


9000 


70 


28 


13 


SFS89, S89 


G330.0+15.0 




1200 


63 


180? 


350 


K96 



a Notcs: Diameters were calculated by using the distances together with the angular 
sizes in Green (2006) catalogue. 
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Fig. 3. — The most likely value of supernova ini- 
tial explosion energy (Eq) derived by compari- 
son with already known distance (d) of the shell- 
type remnants at Sedov-phase equals nearly to 
0.23 x 10 50 ergs. 



Fig. 5. — The most likely value of supernova ini- 
tial explosion energy (Eq) derived by compari- 
son with already known distance (d) of the shell- 
type remnants at Snowplow-phase equals nearly 
to 15.0 x 10 50 ergs. 



1 1 1 1 

E (Most Likely)= 7.0X1 50 ergs 




For Sedov Phase 






(700X10 48 , 0.499) 









1 10 100 1000 10000 100000 

E ( ergs/10" ) 



Fig. 4. — The most likely value of supernova initial 
explosion energy (Eq) derived by comparison with 
already known age (t) of the shell-type remnants 
at Sedov-phase equals nearly to 7.0 x 10 50 ergs. 



E (Most Likely)= 400X1 6 ° ergs 




For Snow-Plow Phase 




" (40000X10*' 


0.868) 







1 10 100 1000 10000 100000 1000000 

E ( ergs/10 4B ) 

Fig. 6. — The most likely value of supernova 
initial explosion energy (Eq) derived by compar- 
ison with already known age (t) of the shell- 
type remnants at Snowplow-phase equals nearly 
to 400 x 10 50 ergs. 
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ison by ages to determine SNe initial kinetic en- 
ergy (Eq) is somewhat meaningful and basically 
no paradox. 

But this measure to determine SNe energy Eq 
through comparison by ages is to some extent 
less reliable than that through comparison by dis- 
tances. Because the initial energy Eq of some in- 
dividual SNRs derived from its age. One can see 
that the obtained remnants distance is more inde- 
pendent to Eq than the age be, and therefore it 
causes less self-contradiction in our work. 

2.2. At radiative phase 

Similarly we list the group of equations for 
shell-type remnants at the third stage (Koyama 
& Meguro-Ku 1987, Kitayama & Yoshida 2005, 
Xu et al. 2005) 

D, c = 1.42( ) 5/2 V (6) 



E(D) - 1.505 x 10~ 19 2 

arcmin 

= 2.88 x 10- u D- c 3 8 n 2 cm - 3 (7) 

V = 10 5 n^/_t( jqJ^ ) 1/8 (8) 

Here, Dp C7 ty T7 n cm -3 7 Siqh z and arcrn i n is 
defined as in section 12.11 as well as their units, 
and tan ( 9arc ™ in ) = ^j^. The fluxes Sighz and 
observational angle 9 arcm i n are already known to 
us for each remnant (table [3]) of which we regard 
these remnants evolving at the snow-plough phase 
since their diameter larger than 36 pc. When the 
initial energy Eq = £ x E^s assumed, then the rem- 
nant diameter D pc (and distance dpc), age t yr and 
electron density be obtained. 

One can see the equations (6) and (8) are rather 
different from equations (1) and (3). But formulae 
(7) and (2) are completely the same. 

For a certain supernova remnant (i) (i = 
1,2,3, ...,20) (table El) and assumed Eq = £ x 
E 48 (£ = 1,2,3, ...,10 5 ) in units of 10 48 ergs, we 
can get the remnant distance d our (^,i)(= d pc ) 
and the age t our (^,i)(= t yr ) by solving the equa- 
tions group above. Then we compare them with 
the already known parameters d true (i) and t true (i) 



(i = 1, 2, 3, 37) listed in table [3] in order to de- 
rive the most likely original energy of SNRs. 

2.2.1. Comparison by distances 

We can compare these resolved distances 
(d OU r(S,,i)) above with the already known ones 
(dtrue(i)) listed in table [3] by 

20 

= ^K m .(^)-<Wz)) 2 

i=l 

n 
i=l 

(e= 1,2,3, ...,10 5 ) (9) 

Figure [5] shows that the most likely value of su- 
pernova initial kinetic explosion energy (Eq) de- 
rived by this method for the S-type remnants at 
Snowplow-phase equals nearly to 15.0 x 10 50 ergs. 

2.2.2. Comparison by ages 

Similarly we can compare these resolved ages 
(iow(£j *)) above with the already known ones 
(ttme(i)) in table [3] by 

20 
i=l 

n 
i=l 

(£ = 1,2,3,...,10 5 ) (10) 

Figure [5] shows that the most likely value of 
supernova initial explosion energy (Eq) derived by 
this way for the shell-type remnants at Snowplow- 
phase equals nearly to 400 x 10 50 ergs. 

2.3. Final results 

From Fig. [3] to Fig. O we get the least value of 
IoqiqEq instead of Eq, therefore one has 

log w E = logioE 02 d + log w E wt 

+log 1 oEQ 3d + log 10 EQ 3 t (11) 

Here, E 02 d, E 02t , E 03d and E 03t corresponding to 
the 4 Eq values in from Fig. [3] to Fig. [BJ Eq is the 
typical explosion energy of shell type remnants. 

Thus we have Eq = 0.99 x 10 51 ergs. 

The publicly accepted value of the SNe initial 
kinetic energy is 1 x 10 51 ergs. 
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3. Discussion and Summary 

To combine theoretical results together with ob- 
servational ones, we have derived the supernova 
average initial explosion energy. Such value is ob- 
viously not the same as that to a individual rem- 
nant. For individual SNR the initial energy is vari- 
ous for different SN explosion events, but it might 
be near an approximation. In Fig. [1] the distri- 
bution of the initial energies of 44 SNRs range 
from 10 49 ergs to 10 52 ergs, mainly concentrate 
to about 10 50 ~ 10 51 ergs. Our 4 theoretical out- 
comes range from 0.23 x 10 50 ergs at minimum to 
400 x 10 50 ergs at maximum. It seems that they 
are to some extent in good consistency regarding 
their divergence. The relative divergence of explo- 
sion energies about j^t is a rather small number 
and acceptable. 

From table [T] one can directly calculate the 
mean initial kinetic energy Eq = 1.7 x 10 51 ergs. 
Or E = 0.85 x 10 51 ergs after excluding SNR 
GO. 0+0.0, namely Sgr A East which owns an ex- 
tremely large initial energy value. These is also 
another method to obtain the initial energy by 
the group of equations combining the physical pa- 
rameters listed in table [2] and table [3) when wc 
regard the age and distance values as the already 
known ones. We gain Eq{i) (i = 1, 2, 3, 37) or 
(i = 1, 2, 3, 20), and the average E is then de- 
rived. Here, in the paper our special method is dif- 
ferent from this, and has provided one more mea- 
sure to derive the SNe typical energy by taking 
somewhat more numbers of the SNRs with inten- 
sified credit. Furthermore, many of the supernova 
remnants in table [2] or table [3] are not the same as 
those in table [TJ Namely the explosion energy Eq 
of rather some remnants in tableland table [3] are 
unknown to us. Because for most of the individ- 
ual supernova remnant in Galaxy their explosion 
energy is rather difficult to derive, here our sta- 
tistical Eq value obtained can be used to compute 
other physical parameters, i.e. the distance, age, 
ISM electron density, magnetic field, etc. of which 
the errors resulted in will be to some extent rather 
small. 

Here we did not distinguish different type su- 
pernova of our adopted sample in the numerical 
analysis. For Type-I and Type-II SN the initial 
explosion energy might be not similar, but one can 
expect that this deviation will be somewhat small. 



The average value of the most likely initial ex- 
plosion energy of supernova remnants thus equals 
to about 0.99 x 10 51 ergs in nice consistence with 
the accepted value. 

JWX thanks J.S. Deng and Y.Z. Ma for their 
assistance and help during the paper work. 
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